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Abstract 

The object of the present investigation is an ensemble of self-avoiding and directed graphs be- 
longing to eight-branching Cayley tree (Bethe lattice) generated by the Fucsian group of a Riemann 
surface of genus two and embedded in the Pincare unit disk. We consider two-parametric lattices 
and calculate the multifractal scaling exponents for the moments of the graph lengths distribution 
as functions of these parameters. We show the results of numerical and statistical computations, 
where the latter are based on a random walk model. 

Keywords: directed random walk; Cayley tree. 

1 Introduction 

A number of stochastic models is taken up to describe the different structures on underlying non- 
trivial geometrical carriers like random environments, percolation clusters, fractals, etc. Restricting 
by physical systems and phenomena, we can note turbulence [TJ [2] , diffusion limited aggregation [3] , 
polymers [H [5], the particle dynamics in a random magnetic field [U [7] as the examples. The last 
decades it seems natural to describe these complicated and chaotic systems by using the multifractality 
concept [U [9] giving us a spectrum of critical exponents which can be connected with the results of 
renormalization group computations (see Ref. |10|). 

We can also find a great number of geometrical objects, chaotic at first sight, within the Lobachevsky 
geometry. Indeed, even limiting by dimension two, there is a possibility to construct the infinite num- 
ber of hyperbolic lattices, existence of which is not allowed in a flat space [11] . This fact inspires us 
in the present paper to turn to Lobachevsky geometry and to demonstrate that the structures, well 
studied in mathematics, might pretend to describe and generate some kind of disorder, like a random 
walk, analytically. Similar idea to consider a directed random walk on a Cayley tree in hyperbolic 
space has been already exploited in Refs. [121113]. We also note that the models of billiards intensively 
use the hyperbolic Riemann surfaces as geometrical carriers [14] 115]. 

In our paper, we shall watch a thread of walker moving along the eight-branching and rooted 
Cayley tree constructed as a dual lattice to the hyperbolic octagonal lattice tiling the Poincare disk 
under the action of Fucsian group. The octagonal lattice cells, whose opposite sides are identified, 
correspond topologically to the two-holed torus, geometry of which is defined here by one independent 
module while other two complex moduli are constrained. (Note that the Riemann surfaces of genus 
g are characterized by 6g — 6 real parameters.) We characterize a walk by spectrum of lengths which 
result from summation of hyperbolic distances between tree sites visited and root point. Introducing a 
partition function and the moments of distribution function, we would like to investigate an ensemble 
of these lengths within the multifractality concept. 

Thus, the model under consideration is deterministic by construction, and partition function is 
actually the truncated series over a symmetry group. To calculate this kind of series, we apply 



both numerical and statistical methods, where the former is needed to exhibit the basic properties of 
the model, which we aim to reproduce by means of the latter approach called for purely analytical 
description. Within the framework of statistical method, we elaborate a scheme of length spectrum 
computation and reduce the partition function to the Markov multiplicative chain in a random walk 
approximation. 



2 Space and Symmetries 

We will deal with a particular case of a two-dimensional Lobachevsky space (with the Gaussian 
curvature K = — 1), namely, a model of the open Poincare disk centered at the origin, 

0= {(x,y) £R 2 \x 2 + y 2 < l} , (1) 

endowed with the metric: 



ds " = 4 m .2_:,2vr ( 2 ) 



dx z + dy l 
(1 — x 2 — y 1 

Mainly, we are interested in the symmetries (isometries) of the metric ([2]), which we will use for 
construction of non-trivial structures on B. However, an application of this metric may not be limited 
by symmetry production only. In Appendix A, we also develop heuristic ideas how to relate the 
curved space to the Euclidean two-dimensional space with potential field. Such a correspondence, 
in our mind, can follow from Maupertuis variational principle, constituting an equivalence between 
particle trajectories (written down below) in these spaces. 

Solution to the geodesic equation in (D,ds 2 ) is the functions: 

cos 4> cosh s + R sin cj> sinh s 

~ Vl + R 2 cosh s + R ' 
R cos 4> sinh s — sin <fi cosh s 

y{s ' = TrTW^shiTR ' 

which are defined in the interval s G (— oo, +oo) and describe an arc inside D with the radius R and 
the center at the point xo = \A + R 2 cos (j), yo = \/l + R? sin0, lying beyond the unit disk. 

Now let us introduce a structure on the Poincare disk. Here we aim to consider an example 
of lattice resulting from operation of a particular Fuchsian group T (stating "periodic law") on the 
(complexified) Poincare disk B [IT]. The associated quotient space B/T is assumed to be a Riemann 
surface of genus g = 2, whose fundamental domain T is hyperbolic octagon in B, which can be viewed 
as the result of gluing two tori together. 

Thus we should first construct the fundamental octagon T and the corresponding Fuchsian group 
r by means of which we tile B by the octagonal lattice. Note that an algorithm connecting J- for 
g = 2 with Fuchsian group T and vice versa is described in Ref. [16J. 

The Riemann surface of genus g = 2 is completely determined by a set of three complex moduli 
corresponding to three corners of the fundamental domain T (in complex plane). The fourth corner 
is then found to adjust the area of the fundamental domain to be 2ir(2g — 2) = 4tt, according to 
the Gauss-Bonnet theorem Applying inversions across the origin to these four corners yield the 
remaining four corners of the octagon. 

For the sake of simplicity, we shall focus on symmetric octagon with the one arbitrary module, 
namely, one pair of the length and angle variables denoted as (a, a). Such an octagon (sketched in 
Fig. 1) can be obtained by deformation of the regular hyperbolic octagon (with a = 2 _1//4 , a = 7r/4), 
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Figure 1: Symmetric octagon with a — 0.8, a — it/3 and generators gk of Fuchsian group. 



well studied in the context of the chaology (see, for example, Refs. [TU [16] and references therein). 
Introducing the complex variable z = x + iy 6 C in order to develop the further analysis in usual 
manner, this construction is as follows. 

Let us assume that the corners of the octagon are at the points aexp (ikir/2), 6expi(a + kir/2), 
where < a, b < 1, k = 0, 3. At this stage, parameter b is unknown and should be found as the function 
of (a, a). Connecting the points by geodesies, we shall require the sum of inner angles (by vertices) of 
the octagon to be equal 2ir in accordance with the Gauss-Bonnet theorem. This requirement allows 
us to determine b and the parameters of geodesies constituting sides of octagon. It is obvious that 
b = a = for the regular octagon (the fixed value of a is required because the homothety property 

is not valid in the space with K ^ 0). 

In the case at a hand, the boundary of octagon, dJ 7 , is formed by geodesies of two kinds (labeled 
by "±" below). These geodesies are completely determined by the radii R± and the angles (j)± + kn/2 
{k = 0,3), defining the positions of the circle centers. Satisfying the conditions imposed above and 
collected in manifest form in Appendix B, we find that 



R± 



1 

2a 



,2\2 



where 



T± = a ± tan ( a 



arctan 



7T 



±1 



(4) 



(5) 



In these formulas it is supposed that < 0+ < a < (f>- < tt/2. Moreover, introducing an angle f3 
by vertices aexp (ikw/2) (the angle by vertices bexpi(a + kw/2) is then equal to ir/2 — /3), 



tan j3 



1 



1 + RiJl + Rl cos; 



1- Jl + RlJl + R 2 _ sin 
(1-a 2 ) 



o 2a 2 cos 2 (a- |) 



2a 2 cos 2 (a -l' 



(6) 



we should control the condition < /3 < tt/2, defining the region of accessibility of parameters (a, a): 

1 



v^cos (a - f ) 



< a < 1, 



(7) 
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which is sketched in Fig. 2. 

To complete the octagon description, the parameter b, pointed out in Fig. 1, can be calculated as 



cos (a — 4> + ) 



1 + Rl 



R\ — tan 2 (a — <f> + ) 



y2a cos (a — 7r/4) 



(8) 



Note that the manifest dependence of the octagon parameters on moduli is necessary in the different 
problems where geometry of the system (carrier geometry) is not fixed. For instance, (a, a) would 
be the (dynamical) variables of evolution parameter in topological gravity describing an universe of 
genus g = 2; it is able to average over (a, a) in statistical physics, etc. 

The Fuchsian group T of a Riemann surface of genus two is generated by four generators go, g±, 
g2, (73 and their inverses. Let generators g^, g^ , acting freely as isometries of Poincare disk, map 
geodesic boundary segments of octagon T onto each other, thereby identifying opposite edges (see 
Fig. 1, right panel). This definition leads to the group relation: 



9o9i 1 929 3 1 9 1 9i9 2 1 93 



id. 



(9) 



Generally, the group of orientation-preserving isometries of (B,(is 2 ) is presented by PSU(1,1) 
SU(1, 1)/{±1} (see textbooks like Ref. p]), where 



u v 
V u 



u, v G C, \u\ 



1 



^(1,1) = 

The group of isometries acts via fractional linear transformations: 

z ^ 7 [ 2 ] = _r_ ; 7 GPSC/(1,1), zGB = {zGC||z| <1}. 



(10) 



(11) 



Then the Fuchsian group T is a discrete group such that T C PSU(1, 1). One also adds the hyperbol- 
icity condition, |Tr j\ > 2 for all 7 G T. Constructing the corresponding Fuchsian group for a given 
octagon, it is convenient to express the generators of T by half turns as follows. 
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Figure 3: The fundamental octagon with a = 0.8, a = tt/3 and the daughter cells. 



Let pk be the mid-point of k-th side, k = 0, 3. Since the opposite sides of octagon have the same 
lengths by construction, the generators are then written as = H{pk) (see Ref. [16]), where 

H(P) = ~t ~TT2 ( 1 + oi P ' 2 i 2 u,2 1 6 PSZ7(1, 1). (12) 



1 — |p| 2 Y 2p 1 + |p| 

Bar over expression, p, means the complex conjugation. 

The operation of matrices H (p) consists in the half turn (rotation with angle it) of geodesic segment 
around the origin z = and the half turn around point p. 

In our model, all pk are parametrized by p± £ D C C: po = p+, Pi = P-, Pi = p+e"^ 2 , and 
P3 = p_e i7r//2 . At this time, p± are functions of two parameters (a, a) and written down below. 
In principal, we could first take pk as "input" and, using then these points, determine the octagon 
parameters. 

To find pk analytically, set 

be ia (l-a 2 ) + ae i7T ^/ 4 (l-b 2 ) . . 

W± = ' (13) 

then one has that 

p ± = (14) 
1 + y 1 — \u±\ 2 

The generators can be directly expressed in the terms of auxiliary variables oj±. Therefore, 
g k = M(u) k ), where 

M(u) = ~^=L= ( l - " 1 € PSU(1, 1), (15) 



\UJ\ 



and wo = w +> oj\ = uj-, 0J2 = w + e 47r / 2 , = w^e" 1 "/ 2 . 

Now it is possible to build the octagonal lattice with a given fundamental octagon. By the con- 
struction, action of the group elements 7 S Y on F produces the daughter cells tiling Poincare disk, 

D = J7J. 

The first eight cells gi^F (k = 0, 3) are shown in Fig. 3. We see that g^FDF C dF as must be. 
Of course, this property is preserved for all neighboring cells (not only for fundamental polygon F). 
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Figure 4: The quarter of fundamental octagon with a = 0.8, a = ir/3, and the Cayley tree sites of first and 
second generations (connected schematically by straight lines but not geodesies). 

3 The Directed Random Walk 

In this Section, we concern with the directed random (self-avoiding) walk (DRW) on a 8-branching 
Cayley tree isometrically embedded in Poincare disk. The tree under consideration is formed by the 
graph connecting the centers of the neighboring octagons (see Fig. 4) and may be physically interpreted 
as a model of directed polymers. 

We start from the general expression for partition function of DRW in disordered environment on 
two-dimensional lattice: 

N 

Zn{q) = ^2X\p{zt-i\zt) exp [qV(z t , t)], (16) 
{z} t=i 

where N is the number of generations; zt is a point of t-th generation of random trajectory; V is a 
random potential of environment; p is a transition probability; parameter q can be mapped into an 
inverse temperature. 

In our model, we consider a deterministic motion on graphs of iV-th generation, sites of which 
correspond to a set of numbers 12, ■ ■ ■ in}, where 1 < it < 8, 1 < t < N. We also assume that 

zt = 7u7 J2 • • • 7i«[0], V(zt, t) = d(0, z t )/N, (17) 

p(z t - 1 \z t )^ P ( it - 1 \i t )=s [ i> , as) 

where {7j|i = 1,8} = {go, gi, g2, 93, 9 \ fi^ 1 ' S^S #3 1 }- B Y construction, points z t and z t -i are the 
centers of neighbouring cells connected by the tree branches. Definition of p(i t -i\it) (consistent with 
definition of generators 73, i = 1,8) excludes a return to the visited sites. Hyperbolic distance d{z,w) 
on (B, ds 2 ) is determined by 

2\z — w\ 2 

COShd(z, W) = 1 + (1 _ N 2 )(1 _ H 2) - ( 19 ) 

In particular case, one has 

d(0, z )=]nii||J. (20) 
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In this way, we come to expression for partition function of our model: 



8 8 / N \ 

Zn(q) = E • • • E ( II 1 ^**- 1 !**) ) exp 

«i=l «jv=1 \t=l / 



1 N 



N 

t 



I2...lt 



(21) 



where Pi lt i 2 ...i t = d(0, 7^7^ ■■■7i t [0]) is the hyperbolic (geodesic) distance between the point of t-th 
generation and the origin O. 

The partition function defined is a sum of Boltzmann weights that depend on an "action integral" 
linear in the hyperbolic distance from a given root point (the origin O), i.e. the length of trajectory: 

N 

L(i 1 ,i 2 , ...,i N ) = ^ft^.-ic (22) 
t=l 

It is worthy to note that Zn(q) is also a function of arbitrary module (a, a), used for obtaining 
the generators of the group T. To explore the properties of a model described by Z^(q\a, a), we shall 
fix the lattice parameters. 

Since the all trajectories are self-avoiding and uniquely determined, the number of graphs Zn(0) is 
equal to the number of end-points, v(N) = 8 x 7 _1 , and the partition function (|2ip can be re-written 
as 

v(N) 

ZsiM) = E eqL ' {N)/N > ( 23 ) 

i=l 

where Li{N) is the integrated hyperbolic length of z-th graph of iV-th generation. 

Thus, the main goal of this work is to evaluate the series (I23h dependent on length spectrum 
{Li(N)}. To resolve this problem, we shall first perform numerical investigations to establish basic 
properties of such a system. Secondly, we shall attempt to calculate the characteristics analytically 
with the use of statistical approach. 

An example of length spectrum obtained numerically is presented in Fig. 5, where the absolute 
value of the number of events depends on particular choice of the width of shell. We see that the 
distribution of lengths Li approaches to the Gaussian distribution, 

(L-Ln) 2 ' 

T 2 



exp(- ^ 2(j ,> I (24) 



with parameters defined by formulas 

v(N) v(N) 

L " = ^v) E L iW, 0% = ^ E (HN) ~ L N )\ (25) 

resulting to Z 5 « 41.97, of 44.5657. 

Although the process of such a convergence is slow, it means that the structure generated by 
the discrete group V satisfactory describes a disordered environment that justifies a random walk 
treatment. Also note that although Fig. 5 was built for the fixed parameters of the lattice, the system 
reveals the same tendency for any admissible pair (a, a). 

We are giving the further analysis of the model within the framework of the concept of multi- 
fractality consisting in a scale dependence of critical exponents [8]. To do this, let us introduce the 
multifractal moments of the order q defined as follows 



,(<?) _ gjyjg) 



< = (26) 
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Figure 5: The number of events as the function of the hyperbolic length L(N = 5) and the Gaussian fit for 
lattice with a — 0.8, a — 7r/3. 



Then the scaling exponent, sometimes called as mass spectrum, is 

where prefactor 2 corresponds to dimension of underlying geometrical carrier (space) and v(N) coin- 
cides with the number of "boxes" tilling structure under consideration. 

Exponent T^(q) is self-averaging and therefore it is reasonable to introduce 

r(q) = lim r N (q), (28) 

N— >oo 

defining the spectrum of fractal dimensions D q = r(g)/(l — q). 

The result of our numerical calculations (limited to N < 5) is shown in Fig. 6. We can conclude 
that the multifractal behavior takes place indeed, that is reflected in non-linear dependence r on q. 

Besides the finding dependence r(q), there is another equivalent way to describe the multifractality. 
It is well known that the set of exponents governing of multifractal moments of type (|26|) is related 
to the spectrum of singularities /(a) of fractal measure. [17j called also the spectral function, which is 
given by the Legendre transform: 

f(a)=qa + r(q), a(q) = (29) 

where a is known as Lipschitz-Holder exponenlQ. In other words, f(a) denotes the dimension of the 
subset characterized by the singularity strength a. In the terms of thermodynamics, a{q) coincides 
with an internal energy e, while /(a) is interpreted as an entropy s(e) [18) . 

The general properties of f(a) are as follows: it is positive on an interval [a m j n , a max ], where 

a min = lim r(q)/(l-q), a max = lim r(q)/(l-q), (30) 

q— »+oo q— > — oo 



1 We preserve conventional symbol for this exponent [8]. It will always be clear from context whether a is an angle 
determining the geometry of octagon or the singularity strength. 
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Figure 6: Multifractal behavior of the exponent r(q) for N = 5. The solid, dashed, and dashed-dotted curves 
correspond to the lattices with (a = 0.8, a = tt/3), (a = 2 -1 / 4 , a — 7r/4), and (a = 0.9, a — 7r/8), respectively. 

and the maximum value of the spectral function gives us the (fractal) dimension of the underlying 
structure. 

Spectral function f(a) obtained for N = 5 is given in Fig. 7. In order to compare curves in Fig. 7, 
we appeal to an information entropy S [8], defined as 

= 1. (31) 

The value of as and f(&s) can De directly calculated at q = 1. 

We find that S(a = 0.9, a = vr/8) ps 1.783 < S(a = 0.8, a = tt/3) rj 1.865 < 5(a = 2- 1 / 4 ,a = 
7r/4) ps 1.95. This relation has a simple treatment: i) due to a special symmetry, the regular octagon 
gives us a highest entropy which cannot exceed 2 in our model; ii) deviation from the regular symmetry 
leads to decreasing the entropy more and more, what says about low probable configurations. 

Since the problem of random walk is mathematically the problem of addition of independent 
random variables, the ultimate behavior of random walkers can be deduced from a central limit 
theorem (CLT). Thus, another approach, applying CLT, can be also developed. 

First, we need to reduce an evolution to the Markov multiplicative process. In order to achieve 
this, let us use the triangle rule in two-dimensional hyperbolic space: 

coshd(0, z t ) = coshd(0, z t -\) coshd(z t , — sinh<i(0, Zt-i) sinhd(zt, zt-i) cosO^t-i, (32) 

where 0t,t-i is angle opposite to side (0, zj). 
Due to identities 

d(0, 7 n [0]) = nd(0, 7 [0]), d(7[0],77[0]) = d(0, 7 [0]), (33) 
which are valid for matrices (fT2|) . (fl~5l) . we get that d(zt,zt~i) = d(Q, 7i t [0]) and 

coshd(0, zt) = cosh d(0, Zt-x) cosh t; lt — sinhd(0, Zt-x) sinh^t cos 0t,t-i, (34) 




o = as = J{ a s), — — 



9 



1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6 2,8 3,0 3,2 

a 



Figure 7: The spectrum of singularities f(a) for N = 5. The solid, the dashed, and the dashed-dotted 
curves correspond to the lattices with (a = 0.8, a = vr/3), (a = 2~ 1//4 ,a = 7r/4), and (a = 0.9, a = 
7r/8), respectively. 

where, accounting for lattice symmetry, = 5f£ + + 

£ + = d(0, 72fc _i[0]), = d(0,72fc[0]), fe = M, (35) 

4 4 

$t = y~!^,2fc-i, #r = y^^,2fc- (36) 
fc=i fc=i 

For relatively large d(0, ^t) and d(0, £t_i), one has 

d(0,z t ) = d(0,z t -i) + ln(cosh£j t — sinh£, t cos 9 t:t -i). (37) 
Now it is useful to introduce 

&j = d(0, 7i7i[0]) " ij -^ ± , i,j = 178, (38) 

describing a logarithmic contribution for the second generation. 

Then, neglecting fluctuations of 9t,t-l i n higher generations, we obtain an estimation for trajectory 
length: 

N 

L(i 1 ,i 2 ,...,i N )^N£ ll +^2(N + l-t)C it _ uit for p(i t -i\i t ) ^ 0. (39) 

t=2 

Note immediately that this formula is exact for i\ = ii = . . . = ijv when #tt_i = 7r. 

In principal, it is also possible to introduce quantities £^\ 2 in ~ d(0, 7^713 • • • 7i„ [0]) playing a r ole 
of higher-order correlations and asymmetric by indices in contrast with symmetric matrix £jj = £j j. 
Moreover, matrices £( n ) are able to account more precisely for values of angles Qt,t-\ or, in another 
word, fluctuations. 

Using Eq. (|39l) . we can easily give the theoretical estimation for the lengths spectrum of iV-th 
generation 

L^ in (N) = Nmm(e+,e-) + N(N ~ ^ minfej), for p(i\j) 0, 
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L ^ M + EE f H l e4^„| ))fe , (4„) 



56 



CI") = iVmax(< + ,<_) + JV(iY 2 '' ^ 



Substituting (|39p in (|2ip . we arrive at the expression for partition function as inhomogeneous 
Markov multiplicative chain: 

8 8 N 

Z N {q) ~ J2--- J2 ^v{qt tl )X\V\i t -i\it), (41) 
ii=l «jv=1 t=2 

where exponent function plays a role of initial condition, and 

N + l-t 



V (i t -i\k) = p(i t -i\it) exp 



N 



(42) 



is transition weight dependent on t. Properties of Eq. (j41j) are discussed in details in Appendix C. 

Thus, we would like to note that the length spectrum {Lj(iV)} and the corresponding partition 
function for any N are completely determined by minimal lengths £± and matrix £jj in a given 
approximation. 

Now, arranging the lengths of graphs in ascending order L± < L% < . . . < L s (s < v(N)), we can 
re-write the partition function ([23]) as 



Z N (q) =v(N)J2 Wie* , £ W* = 1, (43) 



1=1 8=1 



where Zj = Li/N and Wj denotes the weight of the orbit of length Lj. 

Hereafter, coefficients Wj are assumed to be randomly distributed. It means that the value of Wi 
is independent of Wk for i ^ k, allowing the application of the random walk model. 

Next, assuming that the CLT is valid at N — > oo for our system, let us model a distribution of the 
lengths U of iV-th generation graphs by the Gaussian distribution: 

W N (l) = A N ex V (- {l ~f s f y (44) 

where Nl = L N /N, Ns 2 = a 2 N /N 2 when N ->• oo; A N is a normalization constant. We can theoreti- 
cally evaluate I and s 2 by the use of Eqs. ([39]) . ([3D]) . 

We have already seen in Fig. 5 that the Gaussian approximation is in good agreement with nu- 
merical data. So, it is believed that an application of Gaussian distribution may simplify significantly 
the calculations of the series over group T. 

Applying the approximate distribution (144p . we can replace the sums (|43|) with integrals at large 

N: 

<?max(A0 fmax(iV) 

Z N (q) w u(JV) /" VF 7V (/)e <?/ (i/, J W N {l)dl = 1, (45) 

twta(N) imin(N) 
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Figure 8: Multifractal behavior of the exponent r(q) for N = 5 and (a = 0.8, a = 7r/3). The points and the 
curve are obtained on the base of partition functions ([25]) , , respectively. 



where £ m i n (N) = L m i n (N)/N and £ max (N) = L max (N)/N are introduced in order to cut long Gaussian 
tails. 

After simple calculations we obtain that 

Z N {q) « *C N (q) exp N QgV + ql + In 7^ , (46) 

where coefficient 

erf / l max (AQ-M-qiVs 2 \ + erf ( M+qjVs 2 -l min (AQ 



V2Ns 2 J y s/2N^ 



V2Ns 2 J V V2Ns 2 

is responsible for surface effects which are not essential in infinite Bethe lattice; 



X 

2 



erf(x) = -= / e _t dt. (48) 



Comparison of the values of r{q) obtained numerically and calculated on the base of Eq. (|46p is 
presented in Fig 8. We can conclude that a theoretical formula works in the limited interval of q (small 
\q\) and satisfactory describes exponent r{q) for relatively small 1 < iV < oo. However, we can note 
that there is a possibility to improve theoretical result by inclusion of higher-order correlations t^ n \ 
that leads to the processes beyond Markov theory. 



4 Discussion 

Generally, we have concerned here with calculation of the series over Fucsian group as group of 
symmetries of hyperbolic octagonal lattice with one independent complex module. To develop an 
analytical approach, we have considered the model of directed random walk on 8-branching Cayley 
tree as dual lattice for a given octagonal lattice. Then, the truncated series over symmetry group 
was treated as the partition function with Boltzmann weights for ensemble of lengths of tree graphs. 
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Each of these lengths was defined by the sum of hyperbolic distances between the tree sites visited by 
directed walker and the origin. 

As the result, we have analytically derived the spectrum of lengths (|39|) leading to the partition 
function (j41|) in the form of the Markov multiplicative chain allowing us to apply the central limit 
theorem. It says about randomizing environment, where spatial evolution is developed, and applica- 
bility of a random walk method in our calculations. Moreover, we have explored the model within the 
multifractality concept giving us a spectrum of scaling exponents. 

In principal, the investigations performed might be mapped into the problem of calculation of the 
Selberg trace [19] , intensively used for quantization of systems on the Riemann surfaces of genus two 
and requiring the knowledge of geodesies length spectrum (see Ref. [15], for instance). Also, in order to 
obtain the field observables on Riemann surface of the fixed genus and geometry, we can "symmetrize" 
these quantities under the action of corresponding group, that often coincides with the Dirichlet or 
Poincare series computation. We only note a slow convergence of these series, that requires to use the 
accurate approximations which we have touched in this paper. 



Appendix A. Curved Space vs Potential Field 



Let us consider a particle in an attraction potential field U{x,y) > in Euclidean space. The La- 
grangean function and the energy are 



L 



E 



dx\ ( dy\ 



^7 + 



dt J 
dx\ 2 

n) + 



dt J 

dy 
dt 



~-U(x,y). 



Let potential obey the Liouville equation, 

A E lnU = -2KU, 



d 2 d 2 

ox z oy z 



(49) 
(50) 

(51) 



with K = — 1 (negative Gaussian curvature); Ae is the Laplace operator in Euclidean two-dimensional 
space. 

Limiting by the case of central symmetry, when U = U(r), r = y 1 x 2 + y 2 , equation (15ip takes the 
form: 

r/\2 



where U' = dU{r)/dr. 
Solution to Eq. (f52j) is 



U" + -U' 

r 



U(r) 



u 



2U" 



(52) 



(53) 



r 2 sinh 2 (A In r + C) ' 

where A, C are arbitrary constants. 

Following the consequence of Maupertuis variational principle, at the fixed value of the energy E, 
the truncated action function s and the time t are determined from differential relations: 



ds 2 



dt 2 



(2E + U){dr 2 + r 2 dp 2 ), 
dr 2 + r 2 dip 2 
2E + U ' 



(54) 
(55) 
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where (p = arctan (y/x). 

Such a description of dynamical system is based on the equivalence of the trajectories within the 
usual Lagrangean formalism and the geodesies in space with metric (|54|) . 
In the case E = 0, we arrive at the metric: 

ds 2 = U{r)(dr 2 +r 2 d<p 2 ) (56) 

with Gaussian curvature K = — 1. 

Thus we have done transition from the problem of the particle in the potential field in Euclidean 
space to the problem of the free particle in the curved Lobachevsky space. We find that the metric 
([2]) is realized when A = 1, C = 0. 

Since the metric (|56p after replacement r — > exp (— C/A)r, <p —> tp/A gives us ([2]), these metrics 
are isometric. 

Coming back to Eq. (|55p . relation between t (time in Euclidean model) and s (proper time in 
Poincare model now) is derived from expression: 

R 2 Vl + R 2 sinh 




\/l + R 2 coshs + i? LV J 2 J 



1 - R 2 - R) (anli — . (57) 

Then one finds that trajectories for the problem (I49p -(f50l) can be defined in parametric form as 
{*,!!?*(*)} = {t(s),^(s)|s G R}, and 

dx % dx % ( dt\~ 1 



dt ds \ds) ' (58) 



Appendix B. 

In order to describe completely the geometry of fundamental domain T in our model, it is necessary 
to find seven quantities, 

R±, 4>±, b, (3, 7, 

as functions of a and a (see Fig. 1). Parameters R±, <j>± determine the geodesies forming edges of 
octagon, b (or be l ^ a+kn ^ 2 \ k = 0,3) determines the location of four vertices of octagon, j3 and 7 are 
the angles by vertices ae lk7r ^ 2 and be^ a+k7r ^ 2 ^ , respectively. Here we write down seven equations giving 
us the values of these parameters in terms of a and a. 

The first equation results from the Gauss-Bonnet theorem for surfaces of genus two: 

+ 7 = ~- (59) 

Next six equations are derived by using the equations of geodesies in the following (complex) form: 

2 

= R±, k = U73, (60) 

Due to the special symmetry of octagon, it is quite enough to formulate the necessary relations at 
two vertices: z = a and z = be ta . These relations follow from the requirement of intersection of two 
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geodesies, labeled as ±, and are of the form: 



l + a 2 -2aJl + R 2 cos</>+ = 0, (61) 



1 + a 1 - 2a J 1 + Rt sin 0_ = 0, (62) 



a 2 - a ( J I + jRl cos 0+ + J I + i? 2 



sm i 



+J 1 + JK.J1 + Rt sin (0_ - 0+) - R + R_ cos /3 = 0, (63) 



and 



1 + b 2 - 2bJl + R 2 cos(a -</>+) = 0, (64) 



1 + b 2 - 2b^j l + R 2 _ cos - a) = 0, (65) 
6 2 - 6 ^/l + i?2 cos (a - 0+) + ^/l + R 2 _ cos (</>_ - a) J 



+ + + R- sin (0_ - <£+) - cos 7 = 0, (66) 

where (f>+ < a < <p- by construction. 

As it has been noted in main text of the paper, solution to these equations is expressions (3J) — dSJ) 
derived after simple but cumbersome computations. 

Appendix C. 

In this Section, we would like to discuss some properties of the kernel of multiplicative chain (|4ip 
defined at small \q\: 

8 8 N 

K iuiN (q,N) = J2--- E IlP'fe-il**)- ( 67 ) 

12 = 1 «JV — i = l t=2 

First, introducing the vector 

8 

K h (q,N)= E K il>iN (q,N), (68) 

i N =l 

it is worthy to note that this expression allows one the following representation: 

KfaN) = K + (q,N)5f+K_(q,N)67, (69) 

which is obtained on the base of semi analytical observations and reflects the symmetry of the lattice 
in our model; 5^ is defined by Eq. ()36p . Due to this expansion, Eq. (|4ip is re- written as 

Z N (q) ~ 4exp (ql+)K + {q, N) + 4exp (q£-)K_(q, N). (70) 

Next, it is useful to define the probability of transition from initial state i to final state / for N 
steps: 
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In principal, using Pj j(g, N), one can determine a more probable final state / for a given q, N, and 
initial state i. We can also verify that 

P 4i/ (0,iV^oo)^i (72) 

It means that the chain (|4ip is ergodic at q = 0. The same conclusion is obtained by investigation of 
the mean value I (the mean step of graphs) which is introduced in Eq. (|44p and equal to 

r= j&Jt-k (73) 

This result is also independent on initial condition. 
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